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A SPECIAL CASE OF A CONJECTURE BY WIDOM WITH IMPLICATIONS 
TO FERMIONIC ENTANGLEMENT ENTROPY 



ROBERT HELLING, HAJO LESCHKE, AND WOLFGANG SPITZER 



Abstract. We prove a special case of a conjecture in asymptotic analysis by Harold Widom. More 
precisely, we establish the leading and next-to-leading term of a semi-classical expansion of the 
trace of the square of certain integral operators on the Hilbert space L'^(R'^). As already observed 
by Gioev and Klich, this implies that the bi-partite entanglement entropy of the free Fermi gas 
in its ground state grows at least as fast as the surface area of the spatially bounded part times a 
logarithmic enhancement. 



1. Introduction 

In contrast to systems of classical physics, a quantum system composed of two distinguishable 
parts may be in a pure state which is not a product of pure states of its subsystems. Consequently, 
if the total system is in such an entangled pure state, the partial state of each subsystem is not pure, 
in other words, mixed. Following Einstein, Schrodinger, Bell, and others, entanglement may be 
used to rule out interpretations of quantum mechanics which are both local and realistic, similar to 
those of classical (statistical) mechanics. More recently, entanglement has been established as a key 
concept of quantum communication and information theory. For example, quantum teleportation 
and quantum computing heavily rely on it f2','23','3T 40 . 



Partially triggered by the latter theories quantifications of entanglement (e.g., in terms of entropy) 
and consequences thereof are at present intensively discussed for states of many-particle systems. We 
refer to the reviews by Amico, Fazio, Osterloh, and Vedral H^, and by Peschel and Eisler [33]. Here, 
several interesting results and conjectures were put forward. However, in most cases a mathematical 
proof is not yet available and one relies on heuristic arguments, approximate calculations and/or 
numerical observations. This is even true for the entanglement entropy of the ground states of 
quantum spin-chains (see Vidal, Latorre, Ricco, and Kitaev |44]) and of a system being as simple as 
the free Fermi gas. Since in the latter system there is no interaction at all, a non-trivial entanglement 
entropy is solely due to the effective coupling of the particles by the Fermi-Dirac statistics, the 
algebraic statement of Pauli's exclusion principle. 

The interest in entanglement entropy was also sparked from quantum field theory, and in particular 
by toy models for the Bekenstein-Hawking entropy of black holes [sj. Srednicki (SSl and Bombelli, 
Koul, Lee, and Sorkin ^ found numerically that in a semi-classical limit which corresponds to scaling 
the bounded region in d-dimensional Euclidean position space M'' by i? > and taking i? oo, 
the bi-partite entanglement entropy is not a bulk property but scales with the area i?''~"'^|9r2| of the 
boundary surface rather than the volume This so-called area law is thought to be generic for 

field theories with a spectral gap above the ground-state energy. See also the more recent works by 
Cramer, Eisert, and Plenio [17| and by Cramer and Eisert fl6' who proved the area-law scaling for 
harmonic lattice systems. It has been suggested that entanglement might be the mechanism behind 
the black-hole entropy. At first, Bekenstein and Hawking found that black holes behave thermally if 
one interprets the surface gravity as temperature and the area of the horizon as entropy. Especially, 
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there is a "second law" which states that in physical processes the total horizon area can never 
decrease. It is a major challenge for a quantum theory of gravity to show that this is not merely 
an analogy but that the area of the horizon is indeed proportional to physical entropy and to give 
a microscopic explanation thereof. In the framework of string theory this was achieved for extremal 
black holes by Strominger and Vafa 39 and Maldacena and Strominger \29] . More generally, it has 
been argued that the entanglement entropy scales as R'^~^ for d > 2 space dimensions, while for 
d = 1 one expects a logarithmic scaling, Ini?. In a theory with correlation length ^ < oo, heuristic 
arguments suggest that the entanglement entropy stems from correlations across dfl in a layer of 
width ^ > and the absence of long-range correlations is responsible for the area law. However, 
the area law is observed in conformal field theories for 1 + 1 space-time dimensions as well, where 



^ = oo, see Calabrese and Cardy 13 14 



Coming back to simple fermionic systems, Jin and Korepin i25j showed for the first time that 
for free fermions on the one-dimensional lattice Z the entanglement entropy for — [—R, R] OZ 



indeed scales as Ini?, see also Fannes, Haegeman, and Mosonyi 18 . Wolf 50 and later Farkas and 
Zimboras fiol then proved for d> 2 and cubic f2 C Z'' a lower bound on the partial particle-number 
variance that scales as R'^^^hiR. This, in turn, implies that the entanglement entropy grows at 
least as fast as _R''^^ln_R, thereby ruling out an area law. Barthel, Chung and Schollwock [i] and 
independently Li, Ding, Yu, Roscilde, and Haas |28| provided numerical support that the entropy 
itself scales in the same way up to a numerical factor (for d — 2 and d — i). 

To our knowledge, Gioev and Klich 21 were the first to observe an intimate connection between 
the scaling of the entanglement entropy of the free-Fermi-gas ground state and an important conjec- 



ture in asymptotic analysis by Harold Widom 47-49 . This "Widom conjecture" concerns a two-term 



asymptotic expansion of the trace, tr_F(A), for a wide class of analytic functions F of certain integral 
operators A on the Hilbert space L^(]R^), see Equation Q below. The conjecture may be under- 
stood as a multi-dimensional generalization of Szego's asymptotics for Toeplitz determinants and of 
Slepian's spectral asymptotics in classical information theory on the capacity of a communication 



channel which is band limited in both frequency and time. In a similar vein, Gioev 20 established, 
among other things, for the ground state of the free Fermi gas in M'* and rather general C M'' 
with smooth dO, a lower bound on the partial particle-number variance that scales as R'^^^^hiR. 
The main result of the present paper establishes an R'^~^\x\R behavior of that variance itself and 
provides the precise pre-factor in terms of a simple surface integral times a numerical constant. Our 
result is, in fact, a proof of a special case of the Widom conjecture for quadratic F. 

Although we only have a lower bound on the entanglement entropy of the free-Fermi-gas ground 
state, we believe, in accordance with a conjecture by Gioev and Klich ^21j, that this bound refiects 
the correct scaling of the entropy itself up to a numerical factor, which is independent of the Fermi 
sea r C M"* characterizing the ground state, and the region Q, C W^. As they pointed out, their 
conjecture actually goes beyond the Widom conjecture because the entropy corresponds to a non- 
analytic function F (see our Remark |6][iv)). Regardless of the validity of the Gioev-Klich conjecture, 
their works 20 21] were key stimuli to us and apparently also to the authors of (4ll28]. 



The structure of the present paper is as follows: In the next section we formulate the Widom 
conjecture. In Section [3] we prove the Widom conjecture for quadratic polynomials F . Then we 
proceed in Section [4] to compile some background material on fermionic entanglement entropy and 
apply our result. In Section[5]we give an outlook of how to possibly prove the Widom conjecture for 
arbitrary polynomials. The paper ends with appendices on the method of stationary phase, on the 
decay properties of certain Fourier integrals, and on a simple extension of Roccaforte's estimate on 



the volume of certain self-intersections 35 



After having finished the first version of this paper we have learned from Alexander Sobolev 37 



that he has a proof of the Widom conjecture for all polynomials, F, based on pseudo-differential- 
operator calculus. We are grateful for the explanation of his remarkable achievement prior to publi- 
cation. 
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2. The "quadratic" Widom conjecture 

We start with some notation used throughout the paper. If d > 2, we denote a vector t> G M'' 
by a boldface letter, and write v := \v\ := {v ■ v)^/"^ for its norm. Here, we use a dot to denote the 
Euclidean scalar product v w of two vectors w in E**. By A+B := {a+b : a G A,b G B} we denote 
the arithmetic (or Minkowski) sum of a pair of subsets A,BC W^. Wc also write A + b := A + {b} 
for ACM'' translated by 6 G M"^ and RA := {Ra : a e A} for A C M'' multiplied by i? G M. 
For a Borel set A C M'' we denote its volume with respect to the d-dimensional Lebesgue measure 
as |A| := Jf^dx = J^^ dxxAix), where XA stands for the indicator function of A. In particular, if 
A is the positive half-line, Q := X[o,oo [ denotes the right -continuous Heaviside unit-step fimction. 
The Hilbert space of complex- valued, Lebesgue square-integrable functions / : A — >■ C is denoted as 
usual by Li^{K). We use the Bachmann-Landau notation of "little oh" and "big Oh" in asymptotic 
(in)equalities in the sense that for real-valued functions /, g on M, we write 



. f{R)>g{R)+o{h{R)) 
. f{R)=g{R)+o{h{R)) 



if lim n 
if lim/;^oo 
, f{R)^g{R) + 0{h{R)) if linisupj^ 

Next, we formulate our basic assumption. 



f{R)-9(R) > n- 



0; 



f{R)-g(R) _ 
h{R) 

f(R)-g(R) 
h{R) 



< OO. 



ASSUMPTION 1. Letd>2 and C M"* and T C be (C°° -Jsmooth, compact, d-dimensional 
manijolds-with-houndary. The orientation of and of its boundary surface dfl is the one induced 
from R"^ , respectively, from the manifold, and similarly for T and dT . Let a be a smooth, complex- 
valued function on an open set in M'^ x M'^ containing fl x T . 

Note that for such an also the difference 57 — ri:=f2 + (--l)f2 is a smooth, compact manifold- 
with-boundary. For background material in (Riemannian) differential geometry we refer to the 
textbooks [7 10 43 without further notice. 



For two sets 51, F, and a function a as described in Assumption [T] we define for each R > the 
integral operator Aj^ : L^(51) — > L^(f2) by its kernel 



in the sense that 



aR{x,y) (^^) j^dpc^^^^-y^'^ a{x,p) 
[ARf){x):= f dyaR{x,y)f{y), x e n J e L^fl) . 



(1) 



(2) 



Because of 



/ dxidx2 \aB.ixi,X2)\^ < 

Jflxfl 



2d 



\n\^ \r\^ \\a"^ 



oo,si,r 5 



(3) 



where ||a||oo,n,r := sup{|a(a;,p)| 



class, see 



34 



By 



{x,p) G n X r} < oo, the operator Aj^ is in the Hilbert-Schmidt 
Theorem VI. 22(h)], the square Aj^ (and consequently each 



34 



Theorem VI.23] 
natural power A^, fc > 3) is then a trace-class operator. 

We recall that Ajj can be trivially extended to an operator on L^(M'') by viewing _L^(51) as a 
subspace of L^(R'^) and considering xn ^rXXIj where the multiplication operator xn is the orthog- 



onal projection from L 



to L {n). The operators xn^flXn ^-'^d An have the same non-zero 



eigenvalues with the same multiplicities. Therefore, if F is a complex-valued function with F{0) = 
and being analytic on a disc centered at the origin and with radius strictly larger than ||a|loo.f2,r, 
then trF(A«) = trF(xJ^ AflXn). 



4 



R.C. HELLING, H. LESCHKE, AND W.L. SPITZER 



In this and similar situations, Widom 48 , 49 conjectured the beautiful two-term asymptotic 
expansion 

tr F{Ar) ^ [ dxdpF{a{x,p)) (4) 



27r 



2tt 



Ini? / da{x)da{p)\n^ ■ np\F[a{x,p)) + o{R'^ Mni?) 



Here, the linear transformation F ^-^ F is defined by 



47r2 7„ t(l - t) 

Ux E M'' and Up E M."^ denote the outward unit normal vector at cc e dfl, respectively at p e dT, 
and a is the canonical (c? — l)-dimensional area measure on the boundary surfaces dQ and dT. 

Actually, Widom [4^ proved Q in the case that F is a half-space and Q is compact with smooth 
dfl. In 1 35 , Roccaforte considered the case F = M'^ and convolution operators arising from a function 
a not depending on x, and whose Fourier transform is decaying sufficiently fast. Remarkably, he 
proved a three-term asymptotic expansion aR"^ + bR'^~^ + cR'^~^ + o{R'^^^) of tr F{Aji) for certain 
analytic functions F and identified the coefficients a, 6, and c from geometric properties of f2. 

If F{t) = t and if Aji is a trace-class operator then one simply has, 

tr F{Ar) ^ tr Ar ^ ( f dxdpa{x,p)^ f dxdp F{a{x,p)) (6) 

V^^/ Jnxr V^i'/ Jnxr 

for all R > Q. Our main result is the special case of Q with F{t) = as i? — > oo. 

THEOREM 2 ("Quadratic" Widom Conjecture). Under Assumption^ the two-term asymptotic 
expansion 

t^iAnf (7) 

^] I dxdpa{x,pf - -K^ (tt-) Ini? / da{x)da{p) In^ ■ Upl a{x,p)'^ 
27r/ Jnxr 47r^ J Jdnxdr 

+ o{R'^-HnR) 
holds as i? —> oo . 

REMARKS 3. (i) For dimension d ~ I and with fl and F compact intervals, formula ([t]) remains 
true as can be seen by an explicit computation. In this case, the surface integral is simply 
the sum of the four values taken by the function a at the four corners of the rectangle x F. 
(ii) Our proof relies on the method of stationary phase (see |41| ) and an expression for the volume 
of the intersection of a set with its translate as an integral over the boundary (see Theorem 



12 in Appendix |B|) . The proof is elementary in the sense that it does not rely on tools from 
pseudo-differential-operator calculus used by Widom ^49j and recently by Sobolev |37j (see 
the end of the Introduction), 
(iii) In definition ([T]) one could evaluate the "phase-space function" (or "symbol") a more gen- 
erally at {x + \{y — x),p) with A E [0, 1] instead of choosing A = 0. The resulting "A- 



quantization of a" would then lead to an operator Ar x 30 . Here, 1/R plays the role of 
Planck's constant. It can be seen from our proof of Theorem [2J that the asymptotic behavior 
of tr A^ ^ as i? — >■ oo has the same leading term and next-to-leading term as tr A^ for all 
A E [0,1]. 

3. Proof of Theorem [2] 

The proof consists of two parts. The first part deals with the leading term proportional to i?'^ and 
an error term of the order R'^^^ . In the second part we show how the term proportional to R'^^^ In R 
emerges. 
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We start out with a simple change of co-ordinates, u := Xi,v :— Xi — X2, scale v by l/R, and 
hence write 27 p. 524], 



/ dxidx2aR{xi,X2) aR{x2,Xi) 

/ dv duaji{u,u — v) af({u — v,u) 
Jn-n Jn 



d 

dv I dpdqe''"-^P-'i^ 



R 

4^ 



(8) 
(9) 
(10) 



X / dua{u,p)a{u - v/R,q)xn{u) xn{u - v/R) . 



First we expand a{u — v/R, q) at {u, q). The error term is 0{\/R). The integral over u is then of 
the form (e = i?^^) 



/ duf{u)— j duf{u)— / duf{u) 

Jnn(Q+ev) Jn J n\(nn(n+ev)) 



with f{u) := a{u, p)a{u, q). 

Let us define for each a; e the function : M*^ — > C, t> 1— > 7a;(t>) by 

-l^{v) := (271)-'^ / dpa{x,p)e'^-P . 



(11) 



(12) 



Then using the uniform decay, supa,gf2 |7x(''^)| < Cv 2 (see Lemma 11), and Parseval's identity 
we obtain for the "leading term" , 

/ dv ( dpdqa{u,p)a{u,q)e''"<P-i'^ - {2nf [ dpa{u,pf <C/R. (13) 

JBin-n) Jrxr Jr 



For the second term in equation (11 1 we use Theorem 12 with e = R ^. Then, after a change of 
variables we have to analyze the integral 

I ■■= ( ^] I d<j{x) I dv ma.x{0,v ■ n^)'-f^{Rv)j^{-Rv) , (14) 
V^^y Jon Jn-n 

and a remainder term (proportional to w^) which is easy to deal with using the decay of 73,. Namely, 



2d 



dvv^ "/^{Rv)^^{-Rv) 



n-n 



< CR'^-^ 



(15) 



for some constant C. In order to continue with / from ( 14 1, it is convenient to write|^max (0, v ■ n^) = 
Q{v ■ Tlx) v ■ Tlx. Integrating by parts we get 

i2nfv^^iRv) = ^j^dpa{x,p)^e^^^P (16) 

da{p)nj,a{x,p)i}^^ P - [ dp {^a{x,p)) e'^^ A . (17) 

Jr ' 



iR 



For the second integral (over F) we may once more integrate by parts and deduce that it is a term 
of lower order by another factor of R~^. Therefore, for some constant C, 

d 



I + i 



R ^ da{x)da{p) Tlx ■ rip a{x,p) 

Janxdr 

R"^ [ dv e{v ■ Tlx) lA-Rv) e'^^ P 
Jn-n 



(18) 



< CR 



Recall that © is the Heaviside function. 
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The hard part is to analyze the last integral and show that for cr-almost each p e cT one has 



I dvQ{v ■n^)-i^[~Rv)e''^'"P + [2T:i)-'^sgii{n^-np)a{x,p)hiR ^ o{R) 
Jn-n 



(19) 



Here, we need the precise asymptotic expansion ([69| for jx{v) from Lemma 11 and the notation 
employed in its proof. For the resulting phase v i— > R{v-p — v-k) wc are going to apply once more the 
method of stationary phase. To this end, we recall from Lemma [TT] that each k implicitly depends on 
V and use generalized polar co-ordinates, (p, w) , to perform the t>-integration over — fi. In general, 
w G d{^l — fi) and p G [0, 1] are not independent of each other and p does not necessarily cover the 
whole interval [0, 1]. Nevertheless, we may consider the full interval [0, 1] by counting contributions 
with a negative sign if at it) e d{ft — il) the boundary is "inwards" in the sense that w ■ n^, is 
negative. This is sketched in Figure [l] It can be seen that even though pw is not necessarily in 
f2 — SI, points outside are counted with total weight zero while points inside are counted with total 
weight one. 




Figure 1. The integration over 17 — SI. 



Let us pick an ortho-normal basis in M"^ such that the normal vector rip to dV at the given point 
p € dT points in the d-th direction. Locally around p, let dT be given by the graph of a function 
/:[//—> M, /: i—> f{t) with some open Uf C M'^^^; in the notation used in the proof of Lemma 
11 / = /("p'™) and Uf = Un^,m for some m. We write p = {s,f{s)) for some s € Uf and note 
tEat, without loss of generality (by appealing to Sard's Theorem), / is not only critical but has an 
extremum at s. 

In a similar fashion, we can locally write 9(S7 — SI) as the graph of another function g : C/g — > 
M, M I— > g{u) with some open Ug C M*^"^. We assume for the moment that il — fl is convex, or put 
differently that all boundary points are outwards. Then, we partition the integration in ( 19 ) into cones 
V :~ {{pu, pg{u)) : p € [0, 1], u G Ug} C — with some open set Ug C M''"^. Furthermore, instead 
of integrating p over [0, 1] we may integrate over [C / R, 1] for some constant C without changing the 
leading asymptotics of the integral as i? — )■ oo. 

The crucial step is to express the condition (see Lemma[Tof^i)) e := v/v — sgn (v ■ k) rik in these 
new co-ordinates, where, without loss of generality, k = (t^, h{tk)) with h = /(<^'™ ) for some to', 

is now thought of as a fimction of u. Note that = sgn (rip • nk){—^h{tk), 1)/ ^1 + \-§ih{tk)\'^ 
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and thus 



- sgn('U • nfc) sgn [rip ■ rife) 



sgn (5(m)) , = ■ (20) 



— — ■ "to^V" "t>" V-p ■"KJ I 



Let us proceed with the situation that (0, g{0)) G V and call this cone Vq. As there might be several 
disjoint graphs of /("-P'™)^ it is notationally simpler to use their union and call the corresponding 
function again /. This amounts to setting h = f. 

The volume element reads dv = p'^~^f{u) dpdu. As v is parallel to Up = (0, ... ,0, 1) at m = 
and Vq can be chosen small enough, we have Q{v ■ n^,) = 8(sgn(t; • np)np ■ n^), and the only 
contribution to the integral is from those v where g{u) has the same sign as the last component 



of Ux- Using the asymptotics from Lemma 11 we find up to lower-order terms that (using the 
abbreviation /C„ :— IC, , , — r^) 

j dve{v-nxhx{~Rv)e'^'"P (21) 

Jvo 

d + l d+3 f q(u) 
= i(27r) 2 dpp 2 dU ^) 

Jc Ju, (||nP + 3(u)2)(d+i)/4 

X y e((u,.g(n)) .n.)sgn((n,g (n)).n,) ^ 

fet?„ V|det(/,,(t.))| ^ y'^^J^'^))) 

X exp [iRp (m, g{u)) ■ {p - k)] . 



Let us write the last it-integral in the form duip{u) exp [ii?p(/)(M)] . Then we smoothly extend 

to a compactly supported complex- valued function ip on M''"^ and (f) to compactly supported real- 
valued function on M.'^~^ in such a way that^ does not acquire new critical points on the support 
of '0 outside the support of ip- By Proposition [sj this does not change the leading asymptotics of the 
integral. 

Let us investigate now the critical points of the phase function (p : Ug C M''^^ — > M, m 
{u,g{u)) ■ {p — k{u)). Taking derivatives of both sides of equation ([20| and evaluating at it = 
yields (we use the sum convention and sum over indices that appear twice) 



^ = -^^^(^W)^(t.(0))^V (22) 



This implies 



dui 5(0) 
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where f^j ^ denotes the matrix inverse of the Hessian of /. We are now ready to expand the phase 
to second order in it at 0: 

v{p- k)/p (24) 

= u-{s- tk{u)) + g{u){f{s) - f{tk{u))) 

= giO){fis)-fiUiO)) 

1 d^g 



- u. 



-u..,(-^(0) 
dg 



2 duidu 



-(0){f(s) - f(t,m) - i9(0)/4(t.(0))|^|^(0) 



dtrdti 

1 df ^^^^ d'^tr 

2 dtr dujdu-j 



dui du-j 



+ 0{u^). 



Using (231 wc obtain 



V ■ {p-k)/p 

= 5(0)(/(s)-/(tfc(0))) 
+ 0(w'). 



(25) 



u ■ [s 



-tk{0)) 



2 Ui 11 j 



5(0) 



-g.,(0)(/(s)-/(tfc)(0)) 



Then for Vq small enough, the only critical point of the phase function <f> is when tfc(O) = s and 
hence k = p. In this case, 

<l>{u)^p^-^l^^u,u,+0{u^). (26) 
Next, we apply Proposition|9]and conclude that asymptotically (up to next-to-leading terms in 1/ R) 



R'^ dve{v ■ n^)"f^{~Rv)i 

Jvo 



(27) 



i(2^)- 



X exp 



dpp-'gioy 



sgn(n^ • Up) 



V\detiMs))\ 



a{x,p) 



du exp 



' 2g(0) ' 



(27ri) ^sgn{n^-np)a{x,p) / dp p ^ 

JC 

{2ni)^^ sgn{nx ■ Up) a{x,p) \nR. 



Now we address the other cones in the polar decomposition of 17 — and show that they contribute 
to lower order. We choose a parametrization g of the boundary, and cast w e d{fl — Q) in the form 
w = (m, g{u)). Let (f> : d{il — ft) R,w i-^ w-{p — k{w)) be the phase function in the integral ( 19 ) in 
the co-ordinate w. Wc may assume that k{w) ^ p. Since w ■ T^k = 0, a point ioq is a critical point 
of (j) if and only if p— fc(i«o) is parallel to n^jg. In addition, we require that <j){w()) — 0, which means 
that p~ k{wQ) is perpendicular to k;o- For if 0(it)o) 7^ 0, then the /9-integral dp/p exp [ip(j){wQ)] 
would be only of the order 1. For a convex hypersurface d{n — il) these two conditions on Wq cannot 
be satisfied at once and we conclude that only the cone Vq gives the leading contribution. In the 
case of a general hypersurface d{^l — f2), such cases have to be dismissed by the fact they form a set 
of zero surface measure. 

Finally, we discuss the case that — f7 is not convex. By the above argument we only need 
to consider a cone in the direction of Up. Instead of the local parametrization of the boundary 
surface d{fl — Q) as the graph of a single function g we have, in general, a finite number, M , of 
functions 5„ defined on open sets Ug^ for 1 < n < M. We may order them so that Ug^ C C^g„+i 
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and that graph{gM) is the furthest part of the boundary surface, and thus the normal vector points 
outwards. Let us define the cones Vq"'^ := {{pu, pg„ {u)) : p G [0,1], it G Ug^} and let us repeat 
the above calculation for each of these cone. We count their (asymptotic) contributions, namely 
(27ri)~^ sgn(n^ • np)a{x,p)ln R, positive/negative if the normal vector at d{fl — fi) n graph{gn) is 
outwards/inwards. Since is always in the interior of O — i7 and hence M is odd there is only one 
such term that survives this summation and we have finished the proof. □ 

REMARK 4. We have assumed that the boundaries dil and dT are smooth. We believe that our 
proof extends to the case of C^-boundaries. Some regularity, however, is needed as can be seen from 
the example of cubes O = T = [—1,1]'^ with a = 1. In this case the Fourier transform 7, defined in 



(57), is simply given by the product 

,(.)=n^ (28) 

i—1 

with V — {vi, . . . ,Vd). Hence, the leading decay of 7(1') for large \v\ is of the form l^l^" with 

n g {1, . . . ,d} depending on the direction This is in contrast to the leading decay 

in case of a smooth dT. However, the average decay of the Fourier transform is still of the order 



-(''+i)/2 a,s was proved by Brandolini, Hofmann, and losevich 11 for convex sets. In our proof 
of Theorem ^ we critically use the decay behavior for domains F fulfilling our Assumption [l] It is 
obvious that cubes are not covered. However, for the above example of cubes, ([t]) can be proved by 
a direct computation (cf. Remark |3][i)). 

4. Fermionic Entanglement entropy 

We are going to apply Theorem [2] to the ground state of the free Fermi gas in the infinitely 
extended position space R'^. To fix our notation and to supply some background material we first 
consider a slightly more general situation. 

4.1. Entanglement entropy of quasi— free fermionic states. A general system of many fermionic 
particles with separable one-particle Hilbert space H with its scalar product denoted by (•, •), is de- 
scribed by the (smallest) C*-algebra ( [l2j Theorem 5.2.5], [6| Section 4.2]) An generated by the 
unit operator 1 and the annihilation and creation operators a(/) and a*{g) for all f,g G H. These 
operators are bounded and satisfy the usual canonical anti-commutation relations, 

aif)a*ig) + a*{g)aif) = (/,g)l, f,geH, (29) 

a{f)aig) + a{g)a{f) = 0. (30) 

A state p is a linear functional p : A-h C with p{l) — 1 and p{X* X) > for all X G A-h- A 
state p is called quasi-free (and gauge-invariant) [12, p. 43] if there exists a self-adjoint operator D 
on H with < D < 1, such that 

p{a*{f)a{g))^{g,Df), (31) 

and, more generally, 

p(«*(/i)---a*(/™)a(5i)---«(3n)) = { d,t(g„Z?/,.) Ztl (^2) 

for all finite sets {/i, ...,/,„, 171, .. . , (7„} C H. In this sense a quasi-free p is a generalized Gauss- 
ian state, where D plays the role of the covariance. We call D the one-particle density operator 
characterizing p. We note that p is pure if and only if Z) is a projection, that is, = D. 

In order to define the (von Neumann) entropy of a quasi-free state we first introduce the function 

n(t)-l ^ (33) 
'^ '■^\-t\nt-{l-t)\n{l-t) if tG ]0,1[ ■ ^ ^ 

Now, if ri{D) is a trace-class operator then the (von Neumann) entropy, S{p), of the quasi-free p 
characterized by D may be defined as (see [32] Equation (6.9)]), 

S{p) := tr f^{D) . (34) 



10 



R.C. HELLING, H. LESCHKE, AND W.L. SPITZER 



It follows that the entropy of a quasi-free state p is zero if and only if p is pure; this equivalence 
remains true for non quasi-free states but we refrain here from defining the entropy for general states. 

For a general state p and an orthogonal decomposition H = Hi (B H2 into two closed subspaces 
"Hi and 'H2 we use the isomorphism An — Aui ® to define two partial (marginal or reduced) 
states pi and p2 on Aui and Au2i respectively, by 

Pi{X) := p{X®l), X^An,. (35) 

P2{X) := p{l^X), XeAn^- (36) 

Then one has the "triangle" inequality comprised of the Araki-Lieb inequality [s] and the subaddi- 



tivity of entropy 32 Theorem 6.15], 

|5(pi) ~ 5(p2)| < Sip) < Sipi) + S{p2) . (37) 

Here, the left-hand side is zero by definition if S{pi) = S{p2) = 00. As a consequence of ( [37| , the 
partial entropies S{pi) and S{p2) are equal if the (total) state p is pure. A simple quantification 
of the correlations between the subsystems corresponding to T-Li and in the state p of the total 
system, not present in the product state pi® p2, is the (bi-partite) entanglement entropy, 

A5(p) := S{pi ® P2) - S{p) = S{pi) + S{p2) - S{p) > . (38) 

For a pure state p this simplifies to 

AS{p)^2S{pi)^2S(p2). (39) 

In words, for a pure state the entanglement entropy is just twice its partial entropies. 

If p is quasi-free, then pi and p2 are quasi-free, too. More precisely, if p is characterized by D on 
7i as above, then pi {£ E {1, 2}) is characterized by the partial one-particle density operator 

Di EeDEe , (40) 

where Eg : H ^ He is the orthogonal projection from H onto He- Since Di (resp. D2) is the 
zero-operator on H2 (resp. Hi) it is naturally identified with an operator on Hi (resp. H2). By 
construction, the following identities hold. 



pe{a*if)a{g)) = {g,Def), etc. (in analogy to (32 1 ) (41) 



S{pe) = ti r^{De). (42) 

In the special case that the mean of the total number of particles is finite, that is, tv D < 00, then 
the state p is given [12!, Theorem 5.2.14 & pp. 36-37] by a density operator W on the fermionic Fock 
space J-(H) over H. This positive operator of unit trace may be written as 

W = dct(l - D) exp [-^(/„,ln(i?-i - l)/„) a*(/„)a(/m)] , (43) 

n,m 

where {/n} is an arbitrary ortho-normal basis of H. Then one has p2l, 2.5,11, p. 401], [46j 

S{p) = tr T]{D) = -tr W\nW, (44) 



which motivates our definition (34). We stress that trZ) < 00 is not sufficient for tr rj^D) < 00 if 
H has infinite dimension. Conversely, the example D = 1 shows that trri{D) = < 00 is possible 
although trD = 00. 

In the case that (only) tiDi < 00, then (at least) pi uniquely corresponds to a density operator 



Wi on J- (Hi) given by a formula analogous to (43). Accordingly, one then has 



5(pi) = tr77(L>i) = -trT4^ilnW^i . (45) 

Sometimes it is convenient to consider besides the von Neumann entropy also a more general (but 
not subadditive) entropy dating back to Renyi. More precisely, if tr_Di < cxi, we define the partial 
Renyi entropy of order /3 as (cf. [45j Section II. G]) 

5/3(pi) := Intrl^f , /3 € ]0, oo[ \ {1} . (46) 
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Note that S'^(pi) > and lim/^^i Sp{pi) = S{pi). Moreover, the Jensen inequahty hnphes the 
monotonicity, 

(/3-/3')(V(pi)-^Mpi)) >0. (47) 
It may be viewed as a special case of an inequahty between (fractional) absolute moments of a 
random variable, dating back at least to a work of Schlomilch in 1858, see [22j p. 26]. In analogy to 



( 45 ) the quasi-free nature of pi implies 

50(pi)=tr77^(Di), (48) 

where 

m{t):=j^\n{t^ + {l-tf), te[0,l]. (49) 
For later use we also mention the chain of estimates 

2trDi(l - Di) < S2{pi) < (41n2)tr7:>i(l -Di) < S{pi) (50) 
< 5*1/2(^1) < 2trL>y^(l - Di)^/^ < 2trDy^ . 
The first three estimates follow from 2t(l-i) < r]2{t) < (41n2)t(l-i) <?](<) if i G [0,1]. The fourth 



one is (47) with /3 = 1/2 and /3' -> 1, and the last two follow from 771/2(0 < 2t^^'^il - t)^/"^ < 2t^^'^ 
— 1/2 

We now see that trDj^' < cx) is not only sufficient for trDi < 00 but also for tr?7(_Di) < 00. 

While tr Di is physically interpreted as the mean of the number of particles, the quantity tr I?i (1 - 



Di ) occurring in ( 50 1 is the variance of that number in the quasi-free state pi of the subsystem 



corresponding to "Hi. 

4.2. Entanglement entropy of the ground state of the free Fermi gas. Now we consider the 
special case of a free, spinless Fermi gas in d-dimensional Euclidean space M.'^, d € N, at zero absolute 
temperature, that is, in its ground state. In the terminology of Section 4.1 this state p is quasi-free 
and characterized by the Fermi projection D — Q{p — H) on H = L'^{W^). Here, H = h{—i^) is 
a translation-invariant one-particle Hamiltonian given in terms of a smooth "dispersion" function, 
h : M'' M, on momentum space, which tends to infinity near infinity and ensures that _ff is a 
self-adjoint operator on H. The prime example is h{p) = p^, corresponding to the non-relativistic 
kinetic energy (in the absense of a magnetic field). The real parameter p > mi{h{p) : p e M'*} is 
the Fermi energy. Obviously, one has trZ? = 00 but S{p) — tr ri{D) = < 00 due to — D. The 
Fermi sea corresponding to the Fermi projection is given as the lower level set 

F = {p e M'' : h{p) < p} (51) 

in momentum space. 

In order to study the finite-volume properties of the Fermi gas we consider a Borel set 51 C M'' 
with finite volume and thus choose Hi = L'^{il) and T-L2 — i^(M'' \ 57). Then, according to 
Section [4.1[ the partial state pi —: pn of that part of the Fermi gas with bounded position space 
is quasi-free and characterized by 

Di = xh 0(m - H) )ai =■■ Dn . (52) 
We may therefore identify with the operator Ai defined in ^ with the function a — I and F 



given by (51). Moreover, one has (cf. |27, p. 524] for the calculation of tr Dq) 

tvDl < tiDn = (27r)-''|r2||r| < 00 , (53) 



S{pn) = trr]{Dn) < 2trDj/' < °° ■ (54) 



and, by (|50|), even 

Qf \ — T^/' n„ ^ ^ o r^^^ 
1/2 

Here, the finiteness of trD^' and hence that of the partial entropy >5'(pa) of the free Fermi gas in 
its (pure) ground state p, was proved by Gioev and Klich ^21^| by using certain decay properties of 
singular values due to Birman and Solomyak (see also Chang and Ha [l5]). We mention in passing 
that the mean particle density, tr Dn/\n\ = |r|/(27r)'^ is a non-decreasing function of p. 

Theorem [2] has the following 
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COROLLARY 5 (Lower bound on fcrmionic entropy). Suppose T of (51) and VL satisfy Assumption 
[7] Then, the partial entropy S{pj^^) of the free- Fermi- gas ground state satisfies the asymptotic 
inequality 



Sip 



REMARKS 6. 

function a 



> 



In 2 



-V 

2tt J 



Ini? 



da{x)d(7{p) Irix ■ Up] + oiR"^ ^\nR). 



(55) 



dnxdr 



(i) It was already observed by Gioev and Klich 21 that a proof of ([t]) with the 
- 1 would imply (55). As mentioned in the Introduction, Gioev |20l Inequalities 



(1.8) & (1.9)] has previously established a smaller lower bound on S{prq) with the same 
R'^~^ In i?-scaling. 



(ii) An important consequence of the R"^ ^ In i?-scaling of the leading term in (551 is that it rules 
out an area law for the entanglement entropy in the sense that liminfij_>.oo ^'^^d™^^ = oo. 
The /Lt-dependence of that term is encoded in the Fermi surface dT. 

(iii) Gioev and Klich 20,21 also provided an upper bound on S{pun) which is, however, larger 



2, then ( 50 ) 



by an extra factor Ini?. No smaller upper bound is known to us. 
(iv) One may also consider the partial Renyi entropies Sp{pjiQ). For instance, if /3 

gives lower and upper bounds on S2{pRn) in terms of the partial particle-number variance 
tr Duq{1 — -D_Ro), which both scale as R"^^^ In R. More generally, by an informal application 
of the Widom conjecture ^ with a = 1 to the (non-analytic) function F — rj/) from (49) 
and using r^p{\) = (1 + /3)/(24/?) it is tempting to conjecture the exact leading asymptotic 
behavior of the partial Renyi entropy of order /? to be 

1 + /3 f R^ I da{x)da{p)\n^-np\+ o{R'^-^\nR). 



Sf}{pRn) 



(56) 



dnxdr 



24^ 

The von Neumann limit /3 ^ 1 of (56 1 has already been conjectured by Gioev and Klich [21] 
and has stimulated the authors of [4 28 . To our knowledge, the validity of (56) is open 
even for d = 1 and compact intervals and F (cf. Remark |3][i)). See, however, Jin and 
Korepin [25[ Equation (4)] for non-interacting fermions on the one-dimensional lattice Z. 

Proof of Corollary In the (conventional) definition ^ of the operator one keeps Vl fixed and 
(effectively) scales F by R. Here, we need to interchange the roles of the two sets since physically the 
ground state of the Fermi gas in W^, and hence its Fermi sea F is fixed. And one wants to understand 
the asymptotic growth of the entanglement entropy with increasing volume |r2| of the position space 
ri. The required interchangeabilty is justified by the fact that the two products QPQ and PQP in 
terms of two arbitrary orthogonal projection operators Q and P (on L?{W^)) have the same non-zero 
eigenvalues with the same multiplicities. This follows from the singular-value decompositions of QP 
and PQ, see e.g. [36| Section 1.2]. Using the third inequality in ([50| for a lower bound, recalling 
from (53) that ir Aji — trDu^ — (^)'*|ri||F|, and applying Theorem|2 with a = 1 finally gives (55). 

5. Outlook 

Now we show a possible route towards a proof of the Widom conjecture for polynomials of arbitrary 
degree. The reader will have noticed that the essential difficulty is already present for the special 
case a = 1, and that the extension to general a is rather straightforward. In what follows we will 
therefore put a = 1. Then jxiv) of (12) reduces to the simple Fourier integral. 



j{v) := (2^) 



dp e 



iv-p 



It reproduces itself under convolution, that is, 7 * 7 = 7, refiecting the identity Xr 
Proceeding as in equation ([s]) we write for fc e N 

tr (Ar)'' = TT dxj 7(a;j - Xj+i)xRn{xj) , x^+i := Xi , 



Xr- 



(57) 



(58) 
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and introduce new co-ordinates i/q := Xi,yi := X2 — Xi, . . . ,yi^_i 
Vo e e Rfl- yf),. . . , yf._-^ em-y^- y^_^. Then 



:— Xk — Xk~i- Note that 



dyi--- dy^,_^ l{-Vi) ■ ■ ■ l{-Vk-i)l{Vi 
/ dy^ XRn{yo)XRn{ya + Wi) ' ' ' XrUVq + ■ 



■ + yk-i) 

Vk-i) ■ 



(59) 



For the last integral we write (using Lemma 12 ) 

/ dy^ XRn{yo)XRn{Vo + ^i) ' ' ' XRniVo + ■■■ + yu-i) = 

= R'\n\ - R''\n MQnin- yjR) n . . . n (o - yjR - ... - yk-JR))\ 

da{x) max(0, • n^,, . . . , (y^ + . 



R^M-R 



Vk-i) ■ "a 



(60) 
(61) 



+ 0{R -) X XJ^(s2-^^)(yl)■■•Xi^(^2-o)(2/l + •••+yfc-l)• 



For the leading term in (58) we obtain (^)'' |ri||r| by the same argument as for k — 2 with an error 
0{R'^^^); recall that 7 * • • • * 7 = 7. 

Since in ([59]) the function (y^,. . . ,yk-i) ^ li-Vi) ■ ■ ■ l{-Vk-i)l{Vi + • ■ ■ + yfc-i) is symmetric 
we only need to consider the symmetric part of the remaining function, namely of 

(yi,- • ■ ,yfc-i) ^ max(0,yi • n^,, . . . , (yi + . . . + y^,_{) ■ n^) 

X XR(n-U){yi) ■ ■ ■ XR(n~n){yi + ■■■ + yu-i) ■ (62) 

The maximum function by itself can be easily symmetrized by the following quite surprising combi- 
natorial lemma. 

LEMMA 7. Let ai, . . . , a„ he real numbers. Then 

^ max(0, a<^i, a^i + a^2, ■ • • , a^i + • ■ • + a-an) = ^^J max(0, a„i + ... + a„i) , (63) 

where on both sides the summation runs over the n\ permutations of {1, . . . , n} C N. 

The lemma was formulated and used in this version by Widom ^47*, pp. 171,174]. Under the same 
assumptions, Kac 26, Theorem 4.2] presents a proof (due to F. Dyson) that 

n 

max(0, flcri, ao-i + aa2, .■■,a„i + ...+ a„n) = X! X! + • ■ • + «o-fe) , (64) 

o- a k=l 

where Q is (as above) the Heaviside function. It can be easily shown that the right-hand sides of 
(63) and (64) are equal and hence the combinatorial lemma is proved. 

By the transformation kyn we obtain for the power function F{t) — t'^ that 47r^F(l) — — j 



I 



which fits the right-hand side of ( 63 ) . 

Now we come to the next -to-leading term in ( [58| , and consider for 1 < ^ < fc — 1 

c'yi ■ • • c^yfc-i 7(-yi) • ■ ■ 7(-yfc-i)7(yi + • • • + yfc-i) (65) 

X XR(n-Q.{yi) ■ ■ ■ XRin^niVi + ■■■ + yk-i) max(0, (yi + • ■ ■ + y^) ■ n^) ■ 
Then, to leading order, we perform the integration with respect to all variables except v :— yi + 



y^ . This leaves us with the familiar term 



Ir{ 



dv\^{v)\ max(0, t> • ria;) 

'fl(n-n) 

that yields the logarithmic correction term which we know from the k — 2 calculation. 



(66) 
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To complete the proof, one has to show that the error resulting from only symmetrizing the 



maximum function but not the product of indicator functions in ( 62 1 is of lower order as i? — > oo 



Appendix A. Method of Stationary Phase 

We are going to cite two propositions on the method of stationary phase that will be used in 
this paper. To begin with, let us recall that a smooth real-valued function on M.'^~^ has a critical 
point at to € if d(j)(t)/dt\ = 0. Such a point is called non-degenerate if the determinant 

det(j)ij{t) of the Hessian (j)ij{t) := d'^(j){t)/dtitj of is non-zero at t = to- By sgn(j)ij{t) we denote 
the number of strictly positive minus the number of strictly negative eigenvalues of this Hessian at 
t G R''-!. 

PROPOSITION 8. Let r be a smooth complex-valued and let cj) be a smooth real-valued function 
on M"^^^. Moreover, let r have a compact support not containing a critical point of (j). Then 



dtr{t)e'^^^'^ ^0{R-^) (67) 



as R oo for any N Cz N. 



For a proof see (41J Chapter VHI, Section 2, Proposition 4] or 24 Theorem 7.7.1]. The sec- 
ond result (see [4l] Chapter VHI, Section 2, Proposition 6] or [24j Theorem 7.7.5]) deals with the 
asymptotics of the integral in case the phase 4> has a non-degenerate critical point. 

PROPOSITION 9. Let r be a smooth complex-valued and let 4> be smooth real-valued function on 
M.'^^^ . Suppose that cj) has a non-degenerate critical point at to. If r is supported in a sufficiently 
small neighborhood of to (so that there is no other critical points in its support), then there exists a 
sequence {zj)j^fig of complex numbers such that 

e-i«0(to) f ^(i) ei«^(t) ^ i?-(rf-i)/2 I y z^R-^ + 0{R-^) I (68) 

as R^oo for any N eN. In particular, zq = r(to)(27r)('*-i)/2| det0,y(to)p^^^e'T"g"'^'^"(*o). 

In the following K{p) denotes the Gauss-Kronecker curvature of dT at p E dT, and sign(p) 
is the number of strictly positive minus the number of strictly negative eigenvalues of the second 
fundamental form of dT aX p £ dV. 

LEMMA 10. Let T C M'' be a smooth, compact, d-dimensional manifold-with-boundary. Then 
there exists a subset E of the {d — 1) -dimensional unit sphere S'^^^ C R'^ of full Haar measure such 
that for each e £ E there exists a non-empty and finite set ICe C dT such that for all k G /Ce 

(i) e-nfc G {-1,1}, 

(ii) K{k) ^ 0. 

In other words, ICe is the set of non- degenerate critical points of the mapping dT — > M,p i— > p • e. 

As we noted in Remark [i] for the cube T = [—1,1]'', smoothness of dT cannot be relaxed to 
piecewise smoothness without jeopardizing the non-emptiness and/or finiteness of the sets ICe- 

Proof. First we recall the definition of Gauss's spherical mapping G : 9r — S''"^, p i— )■ G{p) := rip, 
and that the curvature K{p) is given by the Jacobian determinant of G evaluated at p G dT. Then 
we define E as E :— G{{p G dT : K{p) ^ 0}). By compactness we have G{dT) = W^^^ and that the 
pre-image ICe ■= G^^({— e, e}) is a non-empty and finite set for each e E E. By Sard's Theorem 
the complement S'^"^ \ is of Haar measure 0. □ 

Choice of co-ordinates: From now on we assume that for fixed e E E the hypersurface dT 
in M'' is locally given by graphs of certain smooth functions : [/g.m ^ defined on some 

open sets Ue,m C M.'^~^ indexed by some to G {1, . . . , M} for some finite M E N {M can be chosen 
independent of e). To be more specific, let us assume for a moment that e — (0, . . . , 0, 1). Then, 
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up to a permutation of the d co-ordinates t = {ti, . . . ,td-i) and Z'-'^''"' we have graph f^^'^'^^ :— 
e R'^-'^ xR : t e Ue,m} and [j^^;^ graph f'-"^"'^ = OF. In addition, we may assume 
that each graph f'-^'"^^ C dV is smaU enough so that it contains at most one k G ICf, (see Lemma 10 ) 
and each fe G /Ce is contained in only one of these graphs. More precisely, for every k G ICe there 
exists a unique rrik G {!,..., M} such that k G graph f'^^'™-''^ and fc = (*fc, some 
tfe G Ue,mk- Note that such a point tk is a critical point of that is, — — = 0. 

t — tfs 

Furthermore, the curvature of 9r at fc G /Ce becomes the determinant of the Hessian of /(•^'™'=) at 
tk (that is, K{k) = det (/^''""^(tfc))), and that sign(fc) = sgn (/^"'"""^ (tfe)). If e = TZ{0, . . . , 0, 1) for 
a suitable rotation 72. then we simply rotate the graphs of /(•^'™) by this TZ. 



LEMMA 11 (Decay of the function 7^). Let T C M'',/Ce C dT, and E C S''"! he as in Lemma\lO 
Finally, let Jxiv) be defined by fTM). Then one has for large v > the asymptotic formula 



jx{v) = -i{2nv)-'i' Y: ''^^||=^«(-,fc)e'-'=+-«)(l + 0(.-^)) (69) 



u/v 

for all v/v G E. The remainder term 0{v^^) is independent of x £ Q. 



Formula ( 69 ) is a slight variant of a standard result that can be found, for example, in ]^A, Theorem 
7.7.14]. For this identification we note that sgn(t; • nfc)cTsig"('=) = ^-f'^ik)^ where "Hormander's 
index" a{k) denotes the number of centers of curvature at fc in the direction v/v minus the number 
of centers of curvature at fc in the direction —v/v. Nevertheless, since we are using these co-ordinates 
in the proof of Theorem [2] we provide a proof based on Propositions [8] and |9] 



Proof. As in (16 1 we use integration by parts to obtain 



V 



IV ^ 



i27rrjUv) = j^- dpaix,p)^e^-P (70) 



A-f / d<7{p)npa{x,p)e''"'P - / dp{^a{x,p))^ 
1''^ ^Jar Jr 



For the second integral in the last equation we may repeat the same integration~by-parts procedure 
with ^ ■ da{x,p) / dp instead of a{x^p). This results in a term of the order v but with the 
same phase as the leading term. Therefore the leading term of 7x(i') as w — > 00 stems from the first 
integral, which we consider in what follows. 

Now, let e := in E be fixed. Moreover, let {tp\)\ be a finite C°°-partition of unity which is sub- 
ordinate to the covering {graph f'-'''"'^)„r of dT in the sense that for each A, supp('0A) C graph f'^^'"^^^''^ 
for some uniquely determined m(A) G {1, . . . , Af}. In addition, if /Cg n supp('0A) is non- 

empty for some A, then we require that ip\{p) = 1 for all p in a neighborhood of this fc(A). 




Figure 2. Co-ordinates for the t-integration. 
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In our co-ordinates we therefore get 

e - f d(j{p)npa{x,p)e''' P = I dtrx{t) e'^*'^*\ (71) 

with certain smooth functions rx and (f>x. In case e — (0, . . . , 0, 1), graphf^'^'^^^ := {{t, Z'^'^'™-' (t)) G 
R'^-^ xR:te Ua.m} as above (see also Figure[2]), and abbreviating p(^)(t) (t, /(^'™(^)) (t)) with 
t G Ue.rn{\)i w6 exphcitly have 

rxit) = sgn(e.np(.,(,))e.(-9/('='™(^»(t)/5t,l)V'A(p('n*))«(^,p'''(*)), (72) 
<Px{t) - e.p(^)(t). (73) 
Formula ( |7T| ) follows from the three facts 
(i) EaV'a(p) = 1, 

(u) dcr{p) = dt^l + |a/(<='™(^))(t)/at|2 for the area measure on graph f^''^"'^^^\ and 
(iii) np = sgn(e • np(A)(t))( - (*)/5*> l)/Vl + |a/('=^™(^))(*)/5t|2 for the unit normal 

vector at p e graph f^'''"^^^^\ 
We note that the signum function in r\ takes either the value 1 or —1 on the whole of Uf,^m(X)_ 



Next, we want to replace the domain of integration Ue.m(\) on the right-hand side of (71) by 
W^^^ without changing the value of the integral. Since ipx has compact support in graph f^^ "^\ we 
smoothly extend rx to M'*"^ simply by setting rx{t) := ii t ^ C^e,m(A)- The phase function (px is 
smoothly extended by Urysohn's Lemma. 



Now, we split the sum in (711 into a sum over those A such that /Ce H supp(V'A) = and those 
for which this intersection is non-empty; in fact, it contains then only a single point, fc(A). Thus we 
have 



V 

V 



d<T(k)nka{x,k)e''''' (74) 

= ^ / dtrA(«)e'^^^W + f dtrx{t)e''"l'^^''\ 

A:K:<,nsupp(V'A)=0 •^^''"^ A:K;ensupp(i/jA)={fc(A)} "'^''-i 

In the first sum we get the decay v^^ for any N according to Proposition 8j In the second sum, 0^ is 



expanded to second order around its critical point tfc(A)- By Proposition 9 we therefore arrive at (69). 
By compactness of i7 we may choose the remainder term 0{v^^) to be independent of a; G fi. □ 

Appendix B. Roccaforte's estimate on the volume of self-intersections 



In 35 Theorem 2.1], Roccaforte proved a theorem which is (by one order of e below) more precise 
than what we need here. See also a previous version by Widom [47| Lemma 2 & 2']). But Roccaforte's 
proof also allows for the inclusion of a smooth function in the integrand. For the convenience of the 
reader we present his proof almost literally and do not claim any originality. Note, however, that 
the derivative of / effects the correction of the order but this is not needed here. 

THEOREM 12 (Roccaforte). Let Q. dW^ he a compact set with C"^ -boundary dfl, Vi, . . . , Vn £ R'^, 
e > 0, and := H (O + Vi) Ci ■ ■ ■ f] {fl + Vn). Let f be a -function defined on fl. Then 

there exists a constant C depending on fl and (the supremum of the derivative of) f so that 



I dxf{x)+e / da{x) f{x) max {OjVk-rix) 

Jn\n,,, ... Jan i<fe<n 



where is the unit outward normal at x Cz dfl. 



<Ce^, (75) 



Proof. Let 17^ :— ^evi,...,evn- Since is compact there exists an Eq and an eg^ubular neighborhood 
A^Eo of such that each x G N^^ can be written uniquely as a; = x — sn^, where x G dfl, is 
the unit outward normal vector at x, and \s\ < Eq. If £ is small enough then 17 \ fi^ C N^^. 
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Let {Uj,ipj} be a finite atlas of co-ordinate neighborhoods covering 9f2. Let Vj := {x G A^^^ : 
ii X = X — srix, then x G Uj}. Define : Vj x M as foUows: if a; = a; — srix G Vj and 

ipjix) = M G M'^^"'^, then 4>j{x) := u + sn, where n := (0, . . . , 0, 1) is the unit vector in W^^^ x M 
normal to K''"^. By the compactness of there exist open sets Nj C Vj such that the Nj are an 
open cover of -/Ve^ and the distance from Nj to the complement of Vj is, for all j, greater than some 
El. If £ is chosen such that max{ jewfe | , 1 < fc < n} < ei then, for all j, k,x G Nj implies x — evk G Vj . 
Let Wj :=Njr]dnc Uj. 

Let {pj} be a partition of unity subordinate to the cover IJ^- Wj of dil. Each extends to a 
function pj on TV, fl O \ by defining pj{x — sn^) :— Pj{x). It now suffices to prove 

/ dx f{x)pj{x) + s / da{x) f{x)pj{x) max (0, ffc • n^) = O(e^) . (76) 

Jn, nn\n^ Jwj i<k<n 

In what follows the index j will be dropped. From the construction of (j) it follows that for y G V, 
y € ilOV ii and only if 4>{y) ■ n = s > 0. Hence for a; G iV, a; G TV H fig if and only if s > and, for 
all 1 < fc < n, (t){x — eVk) ■ n > 0. By Taylor's Theorem a; G n 51^ if and only if s > and, for all 
k, 

(j){x) ■n-e{D^4,){vu)-n + R{e) > , (77) 
where R{e) = O(e^), the estimate being uniform over x since (f> G C'^{V); Dxcf) denotes the derivative 
of at a; with matrix elements {Dx4')ij = '^{x). 

Next we show that for any v G and all a; = S — sn^ G N ^ 

{Dx(j)){v) ■ n = V ■ ng, . (78) 
To see this, let Vt := v — {v ■ nx)nx be the component of v that is tangent to dfl at x. Writing 

{Dx(l>){v) ■ n = {Dx(t>){vt) -n + ivn^) {Dx(j)){n^) ■ n (79) 

it suffices to prove 

{Dx^){vt)-n^O and {Dx^){n^)^n. (80) 

For each fixed sq with |so| < the map <j)^^ : {u, sq) x{u) — Son^iu) describes the hypersurface 
W — SqTI^. The vectors d{x — SQnx)/dui{uQ) form a basis for the tangent space to Vl^ — sqTI^ 
at x{uq) — SQUxiuo). Thus, the derivative D(uo,so)'f'~^ sends vectors {u,0) to vectors tangent to 
W—soTix at x{uq) ~ SQUxiuo) and sends n to n^- Hence it suffices to show that for any x = x — srix, 
the tangent space Tx{W — sn^) — Tx{dfl). But 

d . . d . . s d . , . 
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So, Tx{W — srix) is orthogonal to and hence the same as the tangent space Tx{dVl). 
From ([77| and ([78|, a; G iV n f7 \ 17^ if and only if 

s > and, for some 1 < A: < n, s ~ ev^. ■ + R{e) < . (81) 

Let 

Se ■■= {x e N : s > and for all I < k < n, s - evk ■ ns, + R{e) > 0} , 

and for real 5, 

Is '■— {x N : s > and for all 1 < k < n, s — evk • > 6e^} . 

Then we have 

Nnn\n,^Nnn\s,, (82) 

and the volume of the symmetric difference, 

\S, A Io\ ^ 0{e^) . (83) 



(82) is obvious. To prove (83) note first that since R(e) = O(e^) there is for each (5 > an £5 > so 
that 

e < eg implies Is d Ss d I-s ■ (84) 

Next, 

e < es implies S'^ A /q C I-s \ h ■ (85) 

For, if a; € but x ^ Iq we have by (84) that x G I^s and x ^ lo implies x ^ Is- Similarly, if 

X ^ Ss but X E Iq implies x £ I^s \ Is- From (85 ) it then follows that for e < es, 

|S'e A /o| < 1^-5 \ = / duds \ det D(^a,s)4>~^\ ■ 

J<P{I-s\Is) 

But 

(t){I-s \ Is) ~ {4>{x),x e N : s > 0, —Se^ + e max [v^ ■ n^) < s < Se^ + e max {v^ ■ n^)] ■ (86) 

l<k<n l<k<n 

Thus the above integral is 

< 2fc2 sup I detD(u.s)0"^| / du<25e^M, 

where M := sup{| det ^(q,^)'/'"^! : {u,s) G <j){N)} \(j){W)\ and \(t){W)\ is the {d ~ l)-dimensional 
Lebesgue volume of C M''^^; by the compactness of one can guarantee M < 00. This shows 

This allows us to replace n \ £7^ by iV n \ /q in (76). Changing variables we obtain, 

/ dx{fp)(x)^ ( dMds(/p)o </,-!(«, s) I det7^(^,,)(/.-i I . (87) 

JnnN\io J(p{nnN\io) 

Now we expand (fp) o cj)^^ and | det_D0^^| at s = to first order. Then the last integral equals 

duds [{fp) o 0-i(«, 0) I det i?(a,o)r 'I + 0(5)] , 



where by the definitions of (p and Iq 

(t){nr]N\Io) = {0(a;),a; e N : < s < e max (0, Vfe • n^)} . (89) 

l<A;<n 

Integrating with respect to s yields 

dii(/p)o0-i(M,O)|deti?(^,o)(/'^'| max {Q,Vk-ns,)+0{e^), (90) 

I l<k<n 

which proves our statement by another change of variables. □ 
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